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ABSTRACT 

We present the first measurements of the weak gravitational lensing signal induced by the large 
scale mass distribution from data obtained as part of the ongoing Canada-France-Hawaii Telescope 
Legacy Survey (CFHTLS). The data used in this analysis are from the Wide Synoptic Survey, which 
aims to image ~ 170 square degree in five filters. We have analysed ~ 22 deg 2 (31 pointings) of i' 
data spread over two of the three survey fields. These data are of excellent quality and the results 
bode well for the remainder of the survey: we do not detect a significant 'B'-mode, suggesting that 
residual systematics are negligible at the current level of accuracy. Assuming a Cold Dark Matter 
model and marginalising over the Hubble parameter h £ [0.6,0.8], the source redshift distribution 
and systematics, we constrain as, the amplitude of the matter power spectrum. At a fiducial matter 
density Q m = 0.3 we find as = 0.85 ± 0.06. This estimate is in excellent agreement with previous 
studies. Combination of our results with those from the Deep component of the CFHTLS enables us 
to place a constraint on a constant equation of state for the dark energy, based on cosmic shear data 
alone. We find that wq < —0.8 at 68% confidence. 

Subject headings: cosmology: observations — dark matter — gravitational lensing 



1. INTRODUCTION 

Weak gravitational lensing of distant galaxies by in- 
tervening massive structures in the universe provides us 
with a unique, and unbiased way to study the distribu- 
tion of matter in the universe. Although weak lensing 
has many applications in astronomy, from the study of 
galaxy dark matter halos (e.g., Brainerd et al. 1996; 
Hudson et al. 1998; McKay et al. 2001; Hoekstra et al. 
2004) to galaxy clusters (e.g., Clowe et al. 1998; Dahle 
et al. 2002; Hoekstra et al. 2002a; Cypriano et al. 2004), 
much recent work has been devoted to the measurement 
of the statistical signal induced by large scale structure 
(a.k.a. cosmic shear). 

The reason for the recent popularity of cosmic shear is 
the fact that the signal is a direct measure of the pro- 
jected matter power spectrum over a redshift range deter- 
mined by the lensed sources and over scales ranging from 
the linear to non-linear regime. This straightforward in- 
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terpretation of the signal is rather unique in the growing 
set of tools available for cosmology, and is an important 
feature if we are to determine the values of cosmological 
parameters with high precision. It enables us to mea- 
sure important parameters such as the matter density 
ri m and the amplitude of the power spectrum as, but 
also has the potential to constrain quintessence models 
(Benabed & van Waerbeke 2004; Tereno et al., in prepa- 
ration). The combination of cosmic shear with other well 
understood probes, such as the cosmic microwave back- 
ground (CMB) is particularly powerful (e.g., Contaldi et 
al. 2003; Tereno et al. 2005). 

The measurement of the signal, however, is not with- 
out challenges, which explains why cosmic shear has only 
recently appeared as a useful tool in cosmology. First of 
all, the induced change in the shapes of distant galaxies 
is small (less than 1%), much smaller than the intrinsic 
shapes of the sources themselves. As a result, large areas 
of the sky need to be surveyed in order to reduce the sta- 
tistical errors. The first detections, reported only a few 
years ago (Bacon, Refregier & Ellis 2000; Kaiser, Wilson 
& Luppino 2000; van Waerbeke et al. 2000; Wittman et 
al. 2000), were based on relatively small areas (few deg 2 
at most). Recent developments in the construction of 
wide field imaging cameras on 4m-|- class telescopes have 
made it possible to image much larger portions on the sky 
to warrant accurate measurements of the lensing signal. 
For instance, van Waerbeke, Mellier & Hoekstra (2005) 
presented results based on ~ 12 deg 2 of deep imaging 
data from the VIRMOS-Descart Survey. Other compet- 
itive results are based on shallower data, which cover 
a larger area, such as the Red-sequence Cluster Survey 
(RCS; Hoekstra et al., 2002c) based on 53 deg 2 and the 
75 deg 2 CTIO Lensing Survey (e.g., Jarvis et al. 2003; 
Jarvis et al. 2005). 

The second challenge is the careful removal of observa- 
tional distortions, introduced by the telescope optics and 
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Fig. 1. — Layout of the observations of the Wl and W3 field. Each rectangle corresponds to one of the MegaCam chips. The dashed 
box indicates the planned size of each of the two fields, 72 deg 2 for Wl and 49 deg 2 for W3. Currently 19 deg 2 of Wl and 12 deg 2 of W3 
have been observed. 



atmosphere. Much work has been devoted to deal with 
the point spread function (PSF) and several correction 
schemes have been developed (e.g., Kaiser et al. 1995; 
Bernstein & Jarvis 2002; Refregier 2003) and tested (e.g., 
see Heymans et al. 2005 for the most up to date discus- 
sion). Fortunately the separation of the signal into gradi- 
ent ('E'-mode) and curl ('B'-mode) components provides 
a non-trivial test of the level of residual systematics, in- 
cluding the presence of intrinsic aligments of galaxies or 
observational distortions in the images. In addition, sev- 
eral other tests can be performed to test the accuracy 
of the corrections (e.g., Hoekstra et al. 2002b; Bacon 
et al. 2003; Heymans 2003; van Waerbeke et al. 2005; 
Heymans et al. 2005). The recent results presented in 
van Waerbeke et al. (2005) and Jarvis & Jain (2005) 
have demonstrated that these systematics can be accu- 
rately corrected for, resulting in measurements free of 
'B'-modes. 

The next step is to survey much larger areas on the 
sky to sufficient depth in multiple filters, which enables 
us to probe the evolution of the matter power spec- 
trum. The Canada- France-Hawaii- Telescope Legacy Sur- 
vey (CFHTLS) aims to image ~ 170 square degree in the 
i'-filter down to i' AB = 24.5. This is comparable to the 
depth reached in the VIRMOS-Descart survey, but an or- 
der of magnitude larger in survey area, thus resulting in 
a significant reduction in statistical errors. In addition, 
the fields are observed in 4 additional filters as to ensure 
photometric redshift information for the sources. This 
is an important part of the survey, as it enables tests 
for intrinsic alignments and the study of the evolution 
of the matter power spectrum. The latter significantly 
improves the constraints on cosmological parameters, in 
particular on the equation of state of the dark energy 
(Benabed & van Waerbeke 2004) 

In this paper we present the first results from the 
CFHTLS cosmic shear program, based on 31 pointings, 



resulting in an effective area of ~ 22 deg 2 of i' data, al- 
ready a significant increase compared to previous work. 
The structure of the paper is as follows. In §2 we dis- 
cuss the CFHTLS and the data used in this paper as well 
as the data reduction process. The weak lensing analy- 
sis is described in §3. The resulting cosmic shear signal 
is presented in §4 and the implications for cosmological 
parameters are discussed in §5. 

2. DATA 

The Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) 1 is a joint Canadian-French program to make 
efficient use of Megaprime, the CFHT wide field imager, 
and to address a number of fundamental problems in 
astronomy. Megaprime, equipped with MegaCam, a 36 
CCD mosaic camera with a field of view of ~ 1 deg 2 , 
enables us to obtain deep images of large areas of the 
sky. The survey itself consists of three independent parts, 
each with their own primary science driver. The survey 
has been allocated more than 450 nights over a ~ 5 year 
period. 

The results presented in this paper deal with data col- 
lected as part of the wide synoptic survey, or 'Wide Sur- 
vey' for short. Once completed, it will cover ~ 170 deg 2 
in three patches of 49 to 72 square degrees through the 
whole filter set (u*, g', r', i', z') down to i'=24.5. The 
survey allows the study of the large scale structures and 
matter distribution in the universe through weak lensing 
and the galaxy distribution. The former application is 
the focus of this paper. We also use cosmic shear results 
based on muli-color data from the Deep component of 
the survey, which are described in detail in Semboloni et 
al. (2005). 

For the analysis presented here we use 19 pointings in 
the i'-band from the Wl-field and 12 pointings from the 
W3-field, resulting in an effective survey area 22 deg 2 

1 |http : //www . cf ht . hawaii . edu/ Science/CFHLS 
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after masking. The layout of the current data with re- 
spect to the final survey is presented in Figure ^ These 
data were obtained during the 2003B, 2004 A and 2004B 
observing semesters. In addition, 4 deg 2 of the W2 field 
have been observed, but the area is not contiguous and 
we omit these data from the current analysis. 

2.1. Data Reduction 

Detrended data (de-biased and flatfielded) are pro- 
vided by CFHT to the community through the Cana- 
dian Astronomical Data Centre (CADC). The defend- 
ing is done using the Elixir pipeline developed at CFHT. 
The pipeline also provides photometric zeropoints and in 
most cases a reasonable first order astrometric solution. 

The photometric calibrations are based on observations 
of standard stars during the observing run. These zero- 
points arc only valid under photometric conditions. We 
therefore examine the magnitudes of a large number of 
objects in the images, to check the stability of the pho- 
tometric zero-point and scale the images to the brightest 
image. Most data were taken under photometric condi- 
tions, and the corrections are found to be small. 

The astrometric solutions provided by the Elixir 
pipeline are more reliable for the more recent data, 
whereas a small fraction of chips of the earliest observa- 
tions require a significant revision of the initial astrom- 
etry. The provided astrometric solution, however, is not 
accurate enough for the purpose of our study and needs 
to be improved upon. This process is too time consum- 
ing to be done manually, and is done using an additional 
pipeline. The USNO-A2 catalog could be used to re- 
fine the astrometry, but the number density of sources 
is often too low to warrant stable results. Instead, re- 
trieved a red image from the second generation Digital 
Sky Survey (POSS II) for each pointing. These observa- 
tions have small geometric distortions. The astrometry 
of the POSS II image is calibrated using the USNO-A2 
catalog. SExtractor (Bertin & Arnouts 1996 ) is used 
to generate a catalog of sources with accurate astrome- 
try, with a number density significantly higher than the 
USNO-A2 catalog. In addition, the POSS II images have 
been taken more recently, thus reducing the effects of 
proper motions of the stars. 

This new astrometric catalog is matched to each of 
the MegaCam images. The exposures have been taken 
with different offsets, in order to fill the gaps between 
the chips. We combine the matched catalogs for each ex- 
posure into a master catalog, which contains the average 
positions of the matched objects. This master catalog 
is used to derive the final second order astrometric so- 
lution for each chip. This procedure ensures that in the 
overlapping area, the objects in each exposure are ac- 
curately matched to the same position, which is crucial 
when stacking the images for a weak lensing analysis: er- 
rors in the astrometry lead to additional anisotropies in 
the images. 

Tests show that the resulting astrometric solution for 
each pointing is sufficntly accurate to stack all data into 
a large image. However, in this case different chips can 
contribute to the image at a given position. If the PSF 
properties "jump" between chips, this gives rise to a com- 
plicated PSF anisotropy pattern. Instead, we use only 
those regions of the sky which in the final stacks were 
all observed by the same chip. This avoids complicated 



behaviour of the PSF, at the expense of approximately 
20% of the survey area. In the future we plan to investi- 
gate in more detail to what extent we can deal with the 
PSF on full mosaic images. 

Before stacking the images, we identify stars on chip 22 
(which has the best image quality; see FigureEJ and mea- 
sure the FWHM of the PSF and PSF anisotropy. If the 
seeing of some of the exposures is significantly worse than 
the others, or if the images show extreme PSF anisotropy, 
those images are discarded before stacking the images. 
We stack each chip separately, using the SWarp routine. 
The resulting images are cropped such that the overlap- 
ping regions remain. These images are used in the weak 
lensing analysis presented here. Typically the stacked 
images consist of 7 images, each with an integration time 
of 620s, but because of reasons described above, in a few 
cases we end up with stacks of 6 images. The improve- 
ment in the PSF behavious warrants the minor decrease 
in depth. 

Figure shows the seeing distribution (or PSF size) 
as measured from the stacked images. Panel a shows 
the distribution for stars selected on chip 22, which typ- 
ically has the best image quality. For this particular 
chip we find a median seeing of Of! 71. Towards the edge 
of the field of view the image quality degrades and a 
more accurate representation of the PSF size distribu- 
tion is presented in Figure [SJa, which shows a histogram 
of the distribution based on all chips. In this case the 
median FWHM is 0'.'76. We note that recent changes 
to Megaprime have led to a significant improvement in 
image quality and a large reduction in PSF anisotropy, 
thus reducing the level of systematics in future cosmic 
shear measurements. 

Nevertheless, the image quality in these early data is 
better than 1 arcsecond for all data used here. Figure^ 
shows the PSF size as a function of chip number for the 
best seeing image (filled circles) and worst seeing image 
(open squares). The camera consists of four rows of nine 
chips, which results in the "periodic" changes in the see- 
ing. In both cases the PSF size is best in the middle of 
the camera and increases towards the edges. 

3. WEAK LENSING ANALYSIS 

The weak lensing analysis presented here is based on 
the method proposed by Kaiser et al. (1995) and Lup- 
pino & Kaiser (1997) with a number of modifications that 
are described in Hoekstra et al. (1998) and Hoekstra et 
al. (2000). This particular method has been tested in 
great detail (e.g., Hoekstra et al. 1998, 2002b; Bacon et 
al. 2001; Erben et al. 2001; Heymans et al. 2005) and is 
widely used for cosmic shear studies. Several promising 
alternative methods have been developed recently (e.g, 
Bernstein & Jarvis 2002; Refregier 2003) and tests are 
underway to ensure accurate measurements of the lens- 
ing signal. 

Although the Kaiser et al. (1995) method is not the 
'ultimate' technique to extract the lensing signal from 
the images, it has proven to be one of the most accurate 
techniques currently available when tested on simulated 
data as part of the Shear TEsting Programme (STEP; 
Heymans et al. 2005). In this experiment, various weak 
lensing pipelines were used to measure the shear in sim- 
ulated images. The simulated data consisted of a se- 
ries with 5 different shears (constant across the images) 
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Fig. 2. — (a) FWHM of stars measured from the stacked images of chip 22 for all pointings. This chip typically has the best image 
quality (sec panel c). (b) FWHM measured for all chips, providing a more accurate representation of the PSF size distribution. The arrows 
indicate the median seeing for both distributions, (c) The PSF size as a function of chip number for the best seeing image (filled circles) 
and the worst seeing image (open squares). The camera consists of 4 rows of 9 chips each. The image quality is best in the center of the 
camera, and degrades towards the edge as is indicated by these measurements. 



and 5 different (fairly realistic) PSFs, each introducing 
different systcmatics. Hence STEP provides a test of 
how one can correct for both PSF anisotropy and the 
size of the PSF. However, the simulations do not cap- 
ture all details of real PSFs or imperfections introduced 
by the stacking process. In this experiment, the Hock- 
stra et al. (1998; 2000) implementation of the Kaiser 
et al. (1995) approach was able to recover the lensing 
signal with an accuracy better than 2%. This accuracy 
is sufficient for published cosmic shear results, and ac- 
ceptable for the results presented here, but we note that 
improvements are likely to be needed when analysing the 
complete CFHTLS data set. 

As mentioned earlier, we analyse the stacked images 
for each chip individually. The first step in the analysis 
is the detection of objects, for which we used the hierar- 
chical peak-finding algorithm from Kaiser et al. (1995). 
We select objects that are detected with a significance 
greater than 5<x over the local sky. In addition to the po- 
sitions of the objects, the peak finder also provides fair 
estimates of the object's size. We use this information 
to remove all objects smaller than the PSF. Inspection 
shows that these are either extremely faint objects, or, 
more relevant, spurious detections of diffraction spikes, 
etc. The remaining objects are analysed in more detail, 
which yields estimates for the size, apparent magnitude, 
and shape parameters (polarization and polarizabilitics). 
The apparent magnitudes are corrected for galactic ex- 
tinction using the results from Schlegel et al. (1998). 

The images were inspected by eye to mask out areas 
where the shape measurements could be compromised. 
Potential sources can be cosmetic, such as bleeding stars, 
halos, diffraction spikes, but also astronomical, such as 
HII regions or spiral structure in resolved galaxies. We 
separate the stars and galaxies on the basis of their half- 



light radii. The galaxy shapes are corrected for obser- 
vational distortions as described below. The final, cor- 
rected, catalogs for each chip are combined into a large 
catalog for each field. 

For the weak lensing analysis we use only galaxies 
brighter than i' AB = 24.5, which leaves a sample of 
9.7 x 10 5 galaxies in the Wl field and 6.5 x 10 5 galaxies 
in the W3 field. The surveyed area is ~ 22 deg 2 , result- 
ing in an average galaxy number density of ~ 20 galaxies 
arcmin~ 2 . However, the shape measurement errors are 
larger for faint, small galaxies. Therefore these galaxies 
should be given less weight for an optimal estimate of 
the lensing signal (Hoekstra et al. 2000). 

Instead it is more convenient to define an effective 
number density of galaxies n e fr ; which is related directly 
to the measurement error in the shear for an area of 1 
arcmin 2 , which is the relevant quantity. Hubble Space 
Telescope observations indicate that for well resolved 
galaxies the intrinsic shapes of galaxies results in a dis- 
persion (7 2 ) 1/2 w 0.3 (Hoekstra et al. 2000). If all galax- 
ies were measured 'perfectly', the error in the shear mea- 
surement would be given by a 1 = 0.3/^/n c ff. For the 
CFHTLS data this yields an effective number density of 
~ 12 galaxies arcmin -2 . 

3.1. Correction for PSF Effects 

The observed shapes of the galaxies cannot be used to 
measure the lensing signal, because observational distor- 
tions have significantly altered their shapes in a system- 
atic fashion: PSF anisotropy introduces coherent aug- 
ments in the galaxy shapes and the seeing circularizes 
the images. The key to an accurate measurement of 
the weak lensing signal lies in the adequate correction 
for these systematic effects. As mentioned above, our 
pipeline has been tested extensively (e.g., Hoekstra et al. 
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Fig. 3. — A typical example of the PSF anisotropy as a function of 
position for a stacked Megacam image. The sticks indicate the di- 
rection of the major axis of the PSF and the length is proportional 
to the observed ellipticity of the PSF. In order to show the higher 
order spatial dependence of the anisotropy, we have subtracted the 
average ellipticity. The direction of the average PSF anisotropy is 
indicated in the top right box and the amplitude is indicated as 
well. Although the PSF anisotropy is determined for individual 
chips, the figure clearly shows a large-scale coherent pattern. 



2002b; Heymans et al. 2005) and has been shown to be 
able to recover the weak lensing shear with an accuracy 
of ~ 2%. 

The first step in this procedure is the identification of 
a sample of moderately bright stars which can be used 
to quantify the properties of the PSF (anisotropy and 
'size'). The pattern of PSF anisotropy changes from 
observation to observation, although we found that the 
Megacam PSF is relatively stable, in particular when 
compared to the previous CFH12k camera (e.g., Hoek- 
stra 2004). The PSF anisotropy also varies across the 
field of view and this spatial variation is captured by fit- 
ting a second order polynomial to the shape parameters 
of the stars for each chip. 

We found that a second order polynomial provided an 
excellent fit to the data. Note that the choice of model 
depends on the properties of the camera. For instance 
van Waerbeke et al. (2005) used rational functions (as 
suggested by Hoekstra (2004)) to describe the pattern of 
the CFH12k camera. We intend to improve our charac- 
terization following the procedure developed by Jarvis & 
Jain (2005), which combines information from a series of 
exposures, effectively resulting in a denser sampling of 
the PSF variation. 

A typical example of the PSF anisotropy as a function 
of position across the Megacam mosaic is presented in 
Figure |3] To show the spatial variation in more detail, 
we have substracted the mean anisotropy across the field. 
The resulting pattern is coherent across the field, even 
though the fits were obtained from individual chips. 

Having quantified the PSF anisotropy and its spatial 
variation, we can undo its effect following Kaiser et al. 
(1995) and Hoekstra et al. (1998). The same stars used 



to study the PSF anisotropy are also used to correct for 
the dilluting effect of seeing, as described in Luppino & 
Kaiser (1997) and Hoekstra et al. (1998). 

4. MEASUREMENTS 

To quantify the lensing signal, we measure the ellip- 
ticity (or shear) correlation functions from the galaxy 
shape catalogs. These correlation functions, in turn, can 
be related to the various two-point statistics that are 
commonly used in the literature. The use of the elliptic- 
ity correlation functions also allow for the separation of 
the signal into 'E'-mode (gradient) and 'B'-mode (curl) 
components. Gravitational lensing arises from a grav- 
itational potential and it is therefore expected that the 
lensing signal is curl-free. The amplitude of the 'B'-mode 
then provides a measure of residual systematics. 

In §5 we briefly discuss how these observable two-point 
statistics relate to the matter power spectrum and cos- 
mology. For more detailed discussions we refer to Schnei- 
der et al. (1998) and Bartelmann & Schneider (2001). 
The two ellipticity correlation functions that are mea- 
sured are 



and 



E 



Y,?] WjWj7r,i(Xj) ' 7r,j(Xj) 



(i) 



(2) 



where 9 = |jq — Xj|. 7 t and 7 r are the tangential and 
45 degree rotated shear in the frame defined by the line 
connecting the pair of galaxies. The weights Wi are pro- 
portional to the inverse square of the uncertainty in the 
shear (Hoekstra et al. 2000). For the following, it is more 
useful to consider 



£t(0) = 6t(0)±&r(0), 



(3) 



i.e., the sum and the difference of the two observed corre- 
lation functions. As shown by Crittenden et al. (2002), 
one can derive 'E' and 'B'-mode correlation functions by 
integrating £+(0) and £-(0) with an appropriate window 
function 



where 
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We also present results for other frequently used 
two-point statistics, namely the aperture mass variance 
(M^ p }(9) and the top-hat smoothed variance ("f 2 )(9). 

Of these statistics, the aperture mass is of some par- 
ticular interest, because only for this case are the 'E' and 
'B'-modes uniquely defined, whereas the decompositions 
of the shear correlation function and the top-hat variance 
in E and B modes are defined up to a constant (Critten- 
den et al. 2002; Pen et al. 2002). The aperture mass is 
defined as 
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Fig. 4. — Panel a) 'E' and 'B'-mode shear correlation functions (filled and open points, respectively) measured from ~ 22 deg 2 of i' data 
from the CFHT Legacy Survey. The error bars only indicate the statistical errors determined for the combined signal from the data from 
Wl and W3 fields. Panel b) 'E' and 'B'-mode measurements for the aperture mass statistic. Note that the 'B'-mode is consistent with 
zero on all scales. Panel c) as before but now for the top-hat variance. 



M ap (0) = / dHU(d)K($), 



(6) 



where k is the dimensionless surface density or conver- 
gence. Provided f/(i?) is a compensated filter (i.e., a filter 
such that a constant surface density within the aperture 
yields M ap = 0), the aperture mass can be expressed in 
term of the observable tangential shear 7t using a differ- 
ent filter function Q($) (which is a function of U ('&)), 



M ap (0) = / d^Q^htO?). 



(7) 



We use the filter function suggested by Schneider et al. 
(1998) 
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for 9 < #, and elsewhere. The corresponding Q(9) is 
given by 



6 f-d 2 



#2 



(9) 



for 9 < §, and elsewhere. The 'E' and 'B'-mode aper- 
ture masses are computed from the ellipticity correlation 
functions using 



(M 2 ap )(9) 



and 



(Mi) (9) 



dtid W(%(i?) + W(t?)^_(i>) , (10) 



ddd W(tf)£+(tf)-W(tf , (11 



where W($), and W($) are given in Crittenden et al. 
(2002). Useful analytic expressions were derived by 
Schneider et al. (2002a). Both W(tf), and VV(tf) van- 
ish for > 29, so that (M 2 p ) can be obtained directly 
from the observable ellipticity correlation functions over 
a finite interval. Similarly, the top-hat variance can be 
expressed in terms of the ellipticity correlation functions, 
but with different filter functions, as described in Schnei- 
der et al. (2002a). 

4.1. Lensing signal 

The results for the various two-point statistics are pre- 
sented in Figure The filled points in Figure HJi in- 
dicate the observed 'E'-mode shear correlation function 
as a function of angular scale. The error bars indicate 
only the statistical uncertainty in the measured signal, 
and ignores any contribution from cosmic variance which 
is present in the 'E'-mode signal. However, when esti- 
mating cosmological parameters, as discussed in §5, we 
include estimates for cosmic variance. The open points 
correspond to the 'B'-mode. The signal presented in Fig- 
ure^Jis based on the 22 deg 2 of data for the Wl and W3 
fields. 

As mentioned previously, in order to separate the sig- 
nal into 'E' and 'B'-modes, we have to define a zero-point 
for the 'B'-mode for the shear correlation function and 
top-hat variance. Based on the absence of 'B'-modes in 
the observed aperture mass variance, shown in Figure^, 
the correlated function in panel a has been scaled such 
that the mean 'B'-mode is zero on scales larger than 10 
arcminutes. 

Figure^; shows the top-hat variance (•y 2 ) as a function 
of scale. All three statistics have been used in the past to 
estimate cosmological parameters, and each has its own 
distinct advantage and disadvantage. For instance the 
aperture mass provides an absolute calibration of the 'B'- 
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Fig. 5. — top panel: Aperture mass variance determined from 
Wl (solid points) and W3 (open points) separately. For reference, 
the shaded region indicates the ltr area around the measurement 
from the VIRMOS-Descart survey (van Waerbeke et al. 2005), 
which is of similar depth. The Wl measurements agree well with 
the latter, but the W3 signal is lower than expected. On large scales 
(beyond 10 arcminutes) the agreement between all three measure- 
ments is excellent. Note, however, that the error bars do not in- 
clude cosmic variance, bottom panel: Aperture mass 'B'-mode for 
all three measurements, which are all consistent with no signal. 



mode, but is most sensitive to the matter power spectrum 
on small scales. To probe the largest scales available in 
the surveyed area, we chose to use the top-hat variance 
for our estimates of cosmological parameters, presented 
in §5. 

The results presented in Figure 0] combine the mea- 
surements of the Wl and W3 fields. However, it is use- 
ful to compare the signals obtained from the individual 
fields as well, to check for consistency. The aperture 
mass variances for the two fields are shown in the upper 
panel of Figure El The solid points correspond to Wl, 
whereas the open points are for W3. For reference we 
also show the results from van Waerbeke et al. (2005), 
indicated by the shaded region. The latter uses a similar 
range in apparent magnitude, although the filter is differ- 
ent [Ic vs. i' used here). Hence the amplitude should 
be comparable to our measurements. On scales beyond 
10 arcminutes all three measurements are in good agree- 
ment, but on small scales, the results from W3 appear 
lower than expected. The origin of this discrepancy is 
not clear, and increasing the surveyed area for this field 
might provide a better understanding. The lower panel 
of Figure shows the corresponding 'B'-modes, which 
are all consistent with zero on all scales. 

This paper present the first CFHTLS results based on 
one of two reduction and analysis pipelines. An indepen- 
dent analysis was performed by the Paris group and for 
reference we also show a comparison with their prelimi- 
nary results in Figure [U A detailed discussion of this al- 
ternative analysis will be presented in Fu et al. (in prepa- 
ration) . The top panels in Figure [B] show the aperture 




Fig. 6. — panels a and b: The aperture mass variance from 
the Wl and W3 fields, respectively. The points correspond to the 
results presented here. The shaded regions indicate the la area 
around the measurements from the Paris pipeline (Fu et al., in 
preparation), panels c and d: The aperture mass 'B'-mode signal 
for the Wl and W3 fields. The results from the two pipelines agree 
very well. Note, that the Paris results are based on a smaller survey 
area than the one presented in this paper. 



mass variance for the Wl and W3 fields. The points cor- 
respond to the results from the analysis presented in this 
paper. The shaded regions indicate the la area around 
the measurements from the Paris pipeline (Fu et al., in 
preparation). The lower panels show the corresponding 
measurements of the 'B'-mode signals. The results from 
the two pipelines agree very well. We note, however, that 
the Paris results are based on conservatively masked im- 
ages, resulting in a smaller effective survey area than the 
one presented in this paper. 

4.2. Tests for Residual Systematics 

The absence of a significant 'B'-mode is very encour- 
aging. However, by itself it cannot guarantee that the 
results are free of systematics (the presence of a 'B'-mode 
provides strong evidence for residual systematics, but the 
converse is not true). Fortunately several other tests can 
be performed as well. These tests essentially provide ad- 
ditional checks on the correction for PSF anisotropy, and 
we present results in this section. 

If the correction for the size of the PSF (i.e., seeing) is 
imperfect, this can result in calibration biases (e.g., Hi- 
rata & Seljak 2003). Unfortunately, there is no good test 
of this correction, based on the data alone. Instead, one 
needs to rely on comparison with simulated data sets, 
such as STEP (Heymans et al. 2005). As mentioned 
above, the results of the latter experiment suggests that 
we can recover the shear to better than 2%. We will 
assume a systematic error of this magnitude when esti- 
mating cosmological parameters. 

Depending on the set of cosmological parameters one 
is interested in, one can marginalize over the uncertainty 
in this correction, with little loss in accuracy (Ishak et al. 
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Fig. 7. — top panel: The aperture mass variance computed from 
the uncorrected shapes of the stars. The 'E'-mode (solid points) 
and 'B'-modes (stars) are very similar on scales < 10 arcminutcs, 
but on larger scales the 'E'-mode is larger, bottom panel: The 
aperture mass using the residual shapes of stars after correcting 
for PSF anisotropy. The 'E' and 'B' modes are indistinguishable 
on all scales. The results suggest that the adopted PSF model is 
inadequate on scales between 4 and 20 arcminutes. Note, however, 
that the lensing signal on these scales is ~ 30 times larger than the 
residuals presented here. 



Fig. 8. — The points with error bars show the residual system- 
atics correlation function £ sys as defined in the text. The indicated 
error bars are those of the shear correlation function, and are dis- 
played to show the level of systematics with respect to the ltr 
statistical error. The shaded region corresponds to the ltr region 
around the observed shear correlation itself, which is an order of 
magnitude larger than § sys 



2004). This approach, however, may fail when studying 
the effect of massive neutrinos (e.g., Abazajian & Dodel- 
son 2003). However, such an approach does not take into 
account the variable bias that might arise from pointing 
to pointing seeing variations. Fortunately, as shown in 
Vale et al. (2004), the signal is not affected significantly 
by small variations in calibration errors, thus reducing 
the complexity of the problem for current data sets. 

As was shown in Hoekstra (2004) and van Waerbeke 
et al. (2005) the PSF anisotropy leads to a systematic 
signal which can have different 'E' and 'B'-modes. For 
instance, the systematic 'B'-mode was ~ 3 times lower 
than the systematic 'E'-mode. Furthermore the source of 
an imperfect correction for PSF anisotropy can be two- 
fold. First of all, the correction method itself might in- 
troduce errors. This has been studied in great detail (e.g. 
Hoekstra et al. 1998; Heymans et al. 2005) demonstrat- 
ing that the correction scheme used here, is sufficiently 
accurate. Secondly, as pointed out in Hoekstra (2004), 
an incorrect model for the spatial variation of the PSF 
anisotropy can lead to considerable errors. 

Recently, Jarvis & Jain (2005) have presented an in- 
teresting way to improve the PSF model using the large 
number of exposures taken in a weak lensing survey. 
They analyse the PSF variation using a principal com- 
ponent analysis to identify the dominant PSF anisotropy 
patterns. Under the assumption that only a few param- 
eters are needed to describe most of the variation, this 
approach allows for a much more detailed modelling of 
the PSF. We plan to implement this method for future 
analyses of the CFHTLS data. 

We can test both of these potential points of failure. 
The accuracy of the model describing the spatial vari- 



ation of the PSF anisotropy can be examined by cor- 
recting the measured shapes of the stars and comput- 
ing the aperture mass variance of the residuals. In this 
case, any residual signal is caused by imperfections of 
the model (Hoekstra 2004). Figure0shows the results of 
this test. For reference, the upper panel shows the aper- 
ture mass variance of the stars before correction for PSF 
anisotropy. The 'E'-mode (solid points) and 'B'-modes 
(stars) are very similar on scales < 10 arcminutes, but on 
larger scales the 'E'-mode is significantly larger. Given 
that this signal is larger than the cosmic shear signal 
on scales larger than 10 arcminutes, this figure clearly 
demonstrates the importance of a careful correction for 
PSF anisotropy. 

The bottom panel of Figured shows the aperture mass 
variance using the residual shapes of stars after they 
have been corrected for PSF anisotropy. The 'E' and 
'B' modes are indistinguishable on all scales. The results 
clearly demonstrate that the adopted PSF model is in- 
accurate on scales between 4 and 20 arcminutcs. Note, 
however, that the level of the residuals is about a factor of 
30 smaller than the cosmic shear signal. Although there 
is clearly room for improvement regarding the modelling 
of the spatial variation of the PSF, the results presented 
in Figure [7| suggest that this is not a dominant source of 
error in our current analysis. 

The next step is to examine how well the shapes of the 
background galaxies have been corrected. The amount 
of residual systematics left in the weak lensing signal due 
to imperfect PSF correction can be estimated from the 
correlation between the uncorrected stars and the cor- 
rected galaxies. Bacon et al. (2003) and Heymans et al. 
(2003) defined a useful estimator 
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(e*7) 2 
(e*e*) ' 



(12) 



where e* is the ellipticity of the stars before PSF cor- 
rection and 7 is the shear estimate of the galaxies. The 
estimator is conveniently normalized by the star ellip- 
ticity auto-correlation function, which allows for a direct 
comparison to the lensing signal (j(r)j(9+r)). Note that 
this estimator is sensitive to imperfections in the model 
for PSF anisotropy and imperfections in the correction 
scheme itself. 

The points in Figure |H| correspond to the resulting 
value for £ sys as a function of scale. The indicated error 
bars are those of the shear correlation function, and are 
displayed to show the level of systematics with respect to 
the la statistical error in the lensing signal. The shaded 
region corresponds to the la region around the observed 
shear correlation itself, which is an order of magnitude 
larger than £ sys . These results demonstrate that the 
signal presented in Figure 0] is not significantly affected 
by systematics and can be used to provide reliable con- 
straints on cosmological parameters. Note, however, as 
the CFHTLS progresses the level of residual systematics 
needs to be reduced. 

5. COSMOLOGICAL PARAMETERS 

The observed two-point statistics can be related to the 
matter power spectrum Pg(k), which depends on a range 
of cosmological parameters. For the study presented 
here, the most relevant parameters are the matter den- 
sity Q m and the normalisation erg. In the case of cosmic 
shear, some of the dependence on cosmology also enters 
through the angular diameter distances to the sources. 

5.1. Method 

As discussed above, in order to probe the largest scales 
available in the surveyed area, we choose to use the top- 
hat variance (7 2 ), which can be expressed in terms of the 
power spectrum through 



< 7 2 )(0) = 2tt / dll P K (l) 



Ji(W) 



(13) 



where 9 is the radius of the aperture used to compute 
the variance, and J\ is the first Bessel function of the 
first kind. P K {1) is not the power spectrum itself, but the 
convergence power spectrum, defined as 



™ _ 4c 4 



dw 



(W{w) 
V a{w) 



I 



fK{wY 



(14) 

where w is the radial (comoving) coordinate, wh corre- 
sponds to the horizon, a(w) the cosmic scale factor, and 
Ik(w) the comoving angular diameter distance. W(w) is 
the source-averaged ratio of angular diameter distances 
Dis/Ds for a redshift distribution of sources n(w): 



W(w) = 



dw'n(i 



Jk{w' - w) 



(15) 



f K (w) ■ 

Hence, it is important to know the redshift distribution 
of the sources, in order to relate the observed lensing sig- 
nal to P K (l) and consequently, cosmological parameters. 
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Fig. 9. — Average weight per galaxy as a function of apparent 
'-magnitude. Note that the weights are not normalised 



For our selection of background galaxies (21.5 < i' < 
24.5), the redshift distribution should be very similar to 
the one used by van Waerbeke et al. (2005) in their anal- 
ysis of the VIRMOS-Dcscart data. Currently, only i' im- 
ages have been processed, but as the survey progresses 
photometric redshifts for the sources will be determined, 
which greatly enhances our ability to constrain the cos- 
mology, and to minimize the uncertainty in the source 
redshift distribution. 

Unfortunately, the current observational constraints on 
the redshift distribution are still limited, although a num- 
ber of redshift surveys are targeting the faint galaxies and 
redshift range used for weak lensing (e.g., DEEP2: Davis 
et al. 2003; VVDS: Le Fevre et al., 2004). For the anal- 
ysis here, the best information on the source redshifts 
comes from photometric redshift studies of the Hubble 
Deep Fields (Fernandez-Soto et al., 1999). It is conve- 
nient to parametrize the redshift distribution using 



n(z) = 







l+ct 



exp 



(16) 



The best fit parameters are obtained from a least 
squares fit to the photometric redshift distribution, tak- 
ing into account the adopted weights for the shape mea- 
surements of the sources (see the discussion in §4). The 
weight is predominantly a function of apparent magni- 
tude and the average value is shown in Figure As 
fainter galaxies are, on average, at higher redshift, such 
a weighting scheme will modify the true source redshift 
distribution into an 'effective' one. To compute the best 
fit parameters, we use the Poisson errors in the counts 
as a function of redshift, thus ignoring field-to-field vari- 
ation in the HDF redshift distributions. 

The solid histogram in Figure ^] shows the effec- 
tive source redshift distribution. Comparison with the 
dashed histogram (which corresponds to the unweighted 
case) shows that the weighting scheme slightly lowers the 
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redshift 

Fig. 10. — Effective redshift distribution of galaxies with 
21.5 < Iab < 24.5 determined from the Hubble Deep Field North 
and South. The error bars are the Poisson errors from the finite 
number of galaxies in each bin (note that these do not include 
cosmic variance). The smooth solid curve represents the best fit 
model with a mean source redshift of (z) = 0.81. The dashed his- 
togram corresponds to the redshift distribution when the adopted 
weighting scheme is not taken into account. 

mean source redshift. The smooth curve corresponds to 
the best fit redshift distribution, which has parameters 
a = 1.35, — 1.654 and z s — 0.668. This corresponds 
to a mean source redshift of (z) = 0.81. To quantify 
the uncertainties in the source redshift distribution, we 
keep a and (3 fixed, but vary z s to identify the 68% and 
95% confidence intervals. This yields z s £ [0.632, 0.703] 
with 68% confidence and z s £ [0.613, 0.721] with 95% 
confidence. We marginalise over the latter interval when 
estimating the cosmological parameters. 

The scales probed by the measurement are affected 
significantly by the non-linear growth of structure. As 
shown by Jain & Seljak (1997) and Schneider et al. 
(1998), we cannot use the linear power spectrum, but it 
is necessary to use the non-linear power spectrum. Two 
different approaches to calculate the non-linear power 
spectrum have been proposed and we present results for 
both. The first is based on the scaling formula suggested 
by Hamilton et al. (1991), which was extended to a 
wider range of cosmologies by Peacock & Dodds (1996). 
Peacock & Dodds (1996) provide a prescription to com- 
pute the power spectrum, using a fitting formula which 
is calibrated using numerical simulations. More recently, 
Smith et al. (2003) suggested an approach based on a 
halo model approach to better capture the breakdown of 
the stable clustering assumption in the Peacock & Dodds 
(1996) prescription. We note that both prescriptions are 
based on relatively small numbers of numerical simula- 
tions and that their accuracy is limited. Comparison 
with recent numerical simulations suggest that they are 
accurate to ~ 5%; this number depends on the region 
of parameter space that is probed and larger errors can 
occur (M. White, private communication). 



b 




Fig. 11. — Joint constraints on Q m and ag constraints from 
the CFHTLS Wide data using the Smith et al. (2003) model for 
the non-linear power spectrum. The contours indicate the 68.3%, 
95.4%, and 99.7% confidence limits on two paramaters jointly. We 
marginalised over the Hubble parameter and source redshift distri- 
bution as described in the text. 



5.2. Constraints on O m and as 

Cosmological parameters are estimated by comparing 
the predicted signal m,i to the observed top-hat variance 
di as a function of scale 0i. We consider cold dark matter 
models with a flat geometry (i.e., f2 m + f2A = 1)- We vary 
the parameters of the model, focussing on constraining 
the matter density £l m and the normalisation erg. For 
these parameters we limit the calculations to Q m £ [0,1] 
and a s € [0.5,1.2]. 

As discussed above, we vary the source redshift dis- 
tribution through z s =£ [0.613,0.721], assuming a flat 
prior. Furthermore, the signal depends somewhat on 
the value for the Hubble parameter, for which we use 
h £ [0.6, 0.8] as motivated by the findings of the HST 
Key project (Freedman et al. 2001). In section §5.1 we 
also consider wq, the dark energy equation of state. The 
maximum likelihood function is given by 

C = ( 2 ^)»|C| 1 /2 exp [ (d * " m^idi ~ ™,) T ] • (17) 

Here C _1 is the covariance matrix. C can be decomposed 
as C = C n + C s , where C„ is the statistical noise and 
C s the cosmic variance covariance matrix. The matrix 
C s is computed according to Schneider et al. (2002b), 
assuming an effective survey area of 13.5 deg 2 for the 
CFHTLS Wl and 8.5 deg 2 for W3, a number density of 
galaxies n ga i = 12 per arcmin 2 , and an intrinsic elliptic- 
ity dispersion of a e = 0.3 per component. 

Figure ITT1 shows the joint constraints on n m , a$ using 
the Smith et al. (2003) model for the non-linear power 
spectrum. For reference with other cosmic shear studies, 
we estimate the value of a$ for a fiducial matter density 
of n m = 0.3. For the Peacock & Dodds (1996) model, 
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Fig. 12. — panel a: Dark energy constraints using the measurements from the Wl and W3 fields. The contours indicate the 68.3%, 
95.4%, and 99.7% confidence limits on two paramaters jointly. We marginalised over erg £ [0.7,1.0], h £ [0.6,0.8] and the source redshift 
distribution as described in the text, panel b Results when the measurements from the Deep component (Semboloni et al. 2005) are 
included. We used the Peacock & Dodds (1996) prescription for the non-linear power spectrum. 



we obtain a value of a 8 = 0.88 ± 0.06 (68% confidence). 
The error includes the statistical errors, cosmic variance 
and calibration errors. For the Smith et al. (2003) model 
we obtain a slightly lower value of a 8 = 0.85 ± 0.06. For 
O m < 0.4, the degeneracy between the two parameters 
is well described by as oc fi~ ' 6 . The estimates for cr 8 
based on the Wl and W3 separately also agree well. For 
O m = 0.3, we find a 8 = 0.87 ± 0.07 for Wl and a 8 = 
0.75 ±0.12 for W3. 

These estimates are in excellent agreement with pub- 
lished cosmic shear results from other large surveys. Van 
Waerbeke et al. (2005) list a value of a & = 0.83 ± 0.07 
based on the VIRMOS-Descart survey. Hoekstra et al. 
(2002c) obtained a 8 = 0.861°;^ from the Red-sequence 
Cluster Survey. Jarvis et al. (2005) found a s = 0.8ll";^ 
(95% confidence) from the CTIO lensing survey. The 
constraints obtained here also agree well with those de- 
rived from the CFHTLS Deep survey and a detailed com- 
parison as well as combined constraints with our mea- 
surements are presented in Semboloni et al. (2005). 

5.3. Constraints on dark energy 

One of the major goals of the CFHTLS is the mea- 
surement of the equation of state of dark energy. The 
current lack of (photometric) redshift information seri- 
ously limits the accuracy of such a measurement. Given 
the current limitations of the data, we choose a simple 
model with a constant equation of state: 

P = w p. (18) 

The left panel in Figure IT21 shows the joint constraints 
on Q m and wq based on the measurements of the lens- 
ing signal from the Wl and W3 fields. We marginalised 
over a 8 G [0.7,1.0], h G [0.6,0.8] and the source red- 
shift distribution as described above and considered the 



range of — 2 < < 0. To obtain these results we used 
the Peacock & Dodds fitting formula to obtain the non- 
linear power spectrum. We emphasize that there is not 
yet a reliable analytical fit to the non-linear dark mat- 
ter power spectrum available for non-trivial dark energy 
models. In particular, the halo model proposed in Smith 
et al. (2004) cannot provide an accurate description of 
the wo ^ -1 models because the power spectrum does 
not depend on wq at z — 0. We know that a change in 
wq should affect the background and therefore the change 
in structure clustering, which Smith et al. (2004) do not 
take into account. McDonald et al. (2005) have recently 
extended the halo-model to smaller scales for wo ^ -1 
dark energy models. We will present a detailed dark 
energy measurement in a subsequent work, taking into 
account new model fitting to numerical simulations, and 
the CFHTLS type la supernovae constraints (Van Waer- 
beke et al., in preparation). 

Although we currently lack photometric redshift infor- 
mation, the measurements of the Deep fields presented 
by Semboloni et al. (2005) probe the matter power spec- 
trum at a slightly higher mean redshift, compared to the 
measurements presented here. Hence, we can improve 
our constraints on wq by combining our results to those 
obtained from the Deep survey (Semboloni et al. 2005). 
We stress that this should be considered a consistency 
check rather than a complete dark energy analysis. The 
result of this analysis is presented in the right panel of 
Figure IT21 which gives wo < —0.8 with 68% confidence. 

6. CONCLUSIONS 

We have presented the first cosmic shear analysis based 
on ~ 22 deg 2 (31 pointings) of deep i' imaging data from 
the CFHT Legacy Survey. These observations cover part 
of two of three survey fields and already provide a signif- 
icant increase in area compared to previous work. These 
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early data show a strong variation of the PSF over the 
field of view. We note that recent changes to Megaprime 
have led to a significant reduction in PSF anisotropy. 
Nevertheless, our early results are very encouraging as 
we do not detect a significant 'B'-mode, suggesting that 
the derived lensing signal is free of systematics at the cur- 
rent level of accuracy. Comparison with an independent 
pipeline (Fu et al., in preparation) shows good agree- 
ment and demonstrates that the signal can be recovered 
robustly. 

We assume a cold dark matter model with a flat geom- 
etry and derive joint constraints on the matter density 
f2 m and the normalisation of the matter power spectrum 
as, while marginalising over the Hubble parameter and 
the source redshift distribution. We consider two models 
to calculate the non-linear power spectrum. For a fidu- 
cial matter density of f2 m = 0.3 we find as = 0.88 ± 0.06 
for the Peacock & Dodds (1996) model. Similarly, we 
obtain a value of as = 0.85 ± 0.06 using the Smith et 
al. (2003) approach. These estimates are in excellent 
agreement with previous studies (Hoekstra et al. 2002c; 
Jarvis et al. 2005; van Waerbeke et al. 2005). 

In the coming years we expect to image a total ~ 140 
deg 2 in five filters, allowing us to include photometric 
redshift information for the source galaxies. This will 
greatly enhance our ability to constrain cosmological pa- 



rameters, most notably the equation of state of the dark 
energy (e.g., see forecasts in Tereno et al., 2005). The 
measurements from Semboloni et al. (2005), based on the 
Deep component of the CFHTLS probe a higher mean 
source redshift, thus providing crude redshift leverage. 
When combining our results with the measurements from 
Semboloni et al. (2005) we find w a < -0.8 (68% confi- 
dence) based on cosmic shear measurements alone. 
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